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We propose a forward-backward multiplicity correlation function Cpg, which is experimentally 
accessible, to measure the noncollectivity contribution. It is found that the function is sensitive 
to both of the jet contribution and the small number statistics. We point out that the effect of 
the latter one is also involved in the previous studies of the forward-backward elliptic correlation 
function but was confused as the contribution from jets. We study the CpB in Au+Au collision at 
y/SNN = 200 GeV with a multiphase transport model (AMPT). The result shows that the estimated 
jet contribution is much smaller than previous study due to the finite number statistics which has 
not been noticed before. The connection between this study and the forward-backward elliptic 
correlation function is also discussed. 
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I. INTRODUCTION 

Quantum Chromodynamics predicts a new form of 
matter, a deconfined state of quarks and gluons called 
quark-gluon plasma (QGP). It can be created by heavy 
ion collisions at relativistic energies experimentally 
In a collision, particles are produced by two main pro- 
cesses: soft hadrons arc generated from the hot dense 
medium created by energy sedimentation, while jets are 
produced from initial hard scatterings. Being relevant to 
the two processes respectively, the discovery of the large 
eUiptic anisotropy (^2) [11, SUl and the jet quenching ^ 
have been taken as important evidences of QGP at the 
Relativistic Heavy Ion Colhder (RHIC) 0. 

Particles produced by different processes carry differ- 
ent information of the system. For example, the soft 
hadron azimuthal anisotropy at low transverse momen- 
tum, called elliptic flow, is formed due to the hydrody- 
namic pressure buildup in the initial almond overlap re- 
gion of the colliding nuclei d, Q . At large transverse mo- 
mentum, energetic partons are predicted to lose energy 
by induced gluon radiation Q . This energy loss depends 
strongly on the traversed path length of the propagating 
parton, thus also leads to azimuthal anisotropy 1)2 0, Q . 
The V2 measured from final state particles is the com- 
bination of these two effects. However, the elliptic flow 
V2 and azimuthal anisotropy V2 of jets individually might 
not be the same. Trying to separate the two, previous 
work in this regard has focused on the covariance of 
V2 from forward and backward regions, which has been 
found later that the method is complicated by V2 fluc- 



tuations (lo| . and the jet contribution still can not be 
directly accessed from final state particles. In this paper, 
we approach the problem differently by studying the mul- 
tiplicity correlation. 



II. THE FORWARD-BACKWARD 
MULTIPLICITY CORRELATION FUNCTION 



The procedure is similar to the previous study [9|. To 
insure that jets only contribute to one region while col- 
lective particles contribute to both in a single event, two 
regions with the same size in phase space are chosen. 
Since jets are produced locally, following Ref. 0, we 
make forward-backward rapidity bins with a gap between 
them. The proposed new observable Cp^ is defined as 
the standard correlation function: 
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where Np(^B^ refers to the particle yields in the forward 
(backward) region, and (■ ■ ■ ) means taking the average 
over events. 

To count the contribution of jets to the total multiplic- 
ity, we define a g-parameter: 



(2) 
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Here A^-^(-^) refers to the multiplicity of fiow (jets). Since 
the experimental efficiency docs not bias flow or jet par- 
ticles, forward or backward regions, it is not expected to 
affect the measurement. Similar to 0, in each event, the 
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yields in forward (backward) region can be written as 

NF^{^+rf){N^)+^N^, 

Nb = {1+V^){N^) + {1-0N'', (3) 

with 

N = Nf + Nb = + N'^- (4) 

In the above, N is the total multiplicity in an event, 
and ry^(^) is the event-by-event multiplicity fluctuation 
of flow particles in each region. The jet multiplicity may 
also fluctuates event-by-event, thus no average is taken 
for N'-^ . After taking the average over all events, we will 
get (r;^'^') = 0. For the final state particles, the cor- 
relation between flow particles and jets is believed to be 
small thus ignored, and this is justified by that they are 
produced at different moments during the collision, i.e., 
{N^N'^} = {N^){N'^). £_ describes the jet contribu- 
tion to the forward and backward region. Since jets only 
fall in one side in an event, the assumed value is 1 or 
0. The location is randomly decided, so the probability 
is the same for each side, i.e.,(^) = (1 — ^) = 1/2 while 

Combine Eq. ([2]) and ([3]) and insert into Eq. ([ij, we 
find the relation for g-paramctcr and Cp^: 

C% = -9'+4{rfrj^)il-gr. (5) 

g is the contribution of jets and {rj^rj^) is a measure of 
correlated forward-backward fluctuations. Both of them 
affects the behavior of Cp^ ■ 

Figure [T] shows the relation between Cp^ and the g- 
parametcr. In the ideal case, if flow particles from the 
two regions fluctuate independently, i.e., = 0, as 

assumed in Ref 0, Eq. © can be simply written as 

c% = -g'- (6) 

Then, g can be extracted from the value of Cpg. If the 
system is hydro dominated, 5-^0 and Cpg — 0. While 



if the system is jets dominated, <? — 1 and Cpg — > — 1. 
On the other hand in a realistic case, if flow parti- 
cles from the two regions do not fluctuate independently 
{{rj^rj^) 0), when g is small, Cp^ will be larger than 0, 
and when g is large, the effect of (Tj^Tj^) can be ignored. 

The Eq. ([5|) works in the low px with relatively large 
number statistic. However, if there is only one flow parti- 
cle produced a particular px bin in a event, it can only fall 
into one side and behaves like a jet. In that case, the mul- 
tiplicity correlation function will automatically decrease 
when the particle yield is much smaller than 1. This ef- 
fect has not been discussed in previous studies d, [lo| . 
and can be easily confused with the contribution from 
jet. To explore this effect, a simulation was performed 
as followed. The flow particles were generated with a 
exponentially decreases distribution, and randomly fell 
into forward/backward region. Adopting the consistent 
terminology as in Q , we call them " flow-only" particles 
throughout this paper. The red curve in Fig. [5] shows 
the result of flow-only simulation. It can be seen that 
although there is no jet, Cpg is still approaching to —1 
at large pt- This phenomena is solely due to, for a given 
Pt bin, the exclusion of events for which both Np and 
Nb for that pt bin are zero, from the CpB calculation. 
This happens even when the two samples are completely 
independent of each other. 

With the simplest assumption that when particle yield 
N{j)t) is smaller than 1, maximum one particle can be 
found in either forward or backward region, the proba- 
bility of finding an empty pT bin is (1 — N{pt)/2)'^ ■ If 
one calculates ^jv^^ffi^^ without the exclusion of events 
with empty bins as mentioned above, for two independent 
samples it should be at unity. Considering the exclusion, 
one needs to scale each averaging term in by 

a factor of ]^_(x_jv'(p^)/2)^ ' then Cpg^px) can be written 
as: 

C^BiPT) = -il-NipT)/2r. (7) 
Thus even for independent samples, with the exclusion 
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of empty bins, Cpg{pT) is decreasing with increasing 
N{pt). The black dash hne in Fig. His a plot of Eq. ©. 
Wc can sec that when pT is larger than 2.2 GeV/c, the 
equation can well describe the simulation result. While 
when pt is smaller than 2.2 GeV/c, the simulation result 
is lower than the black dash curve. This is because when 
particle yield is not so small, there might be more than 
one particle generated in the forward/backward region 
for the flow-only simulation. Therefore, from the Eq. (|6]) 
and Fig. [5] we can conclude that both of the jets and 
the scarcity of particles can cause the decreasing trend of 
Cpg{pT)- We may also expect that it will decline even 
faster when jets are involved in the particle-rare case. 
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III. AMPT MODEL 

The forward-backward multiplicity correlation func- 
tion Cpg is studied in a multiphase transport model 
(AMPT) There are four main components in 

this transport model: the initial conditions, the parton- 
parton interactions, the conversion from the partonic 
to the hadronic matter and the late hadronic inter- 
actions. The initial conditions are based on the HI- 
JING model [l^ in which the eikonized parton model 
is employed. It includes the spatial and momentum 
distributions of minijet partons from hard processes 
and strings from soft processes. The parton-parton in- 
teractions and the time evolution of partons is then 
treated by the Zhang's Parton Cascade (ZPC) [l^ model. 
The hadronization process is described by a combined 
coalescence and string fragmentation model. A rel- 
ativistic transport (ART) model which includes 

baryon-baryon, baryon-meson and meson-meson elastic 
and inelastic scattering is employed to describe the late 
hadronic process. In our study, we analyzed the events 
from string melting AMPT model for Au -I- Au collisions 
at ^/snn = 200 GeV with parton cross sections equal to 
10 mb. Pseudorapidity regions [—1,-0.5] and [0.5, 1] are 
chosen as forward and backward region respectively. 

The Pt dependence of multiplicity correlation function 
is shown in Fig. [31 We can see that Cpg almost stays as 
a constant when px < 1-5 GeV/c. Since particles from 
high energy jets give only a small contribution at those 
Pt, the value of Cpg reflects the flow multiplicity fluc- 
tuation, A{r]^r]^). Although the particle yield decreases 
exponentially, A{rj^rj^) shows little pT dependence for 
this region. To obtain the jet contribution, Cpg is fitted 
by the function: 



C%^^q{pTf+^{v''if){l-q{PT)) 
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and q{pT) is parameterized using a tanh function to de- 
scribe threshold behavior between the two regimes, as 
suggested in [§]: 



q{PT) = {1 + tanh[(pT - Pc)/Pw]}l2. 



(9) 



Here q{pT) is not simply g{pT) anymore. Instead, it is 
a convolution of both g{pT) (dominated at low pt) and 



FIG. 3: (Color online) Cpg as a function of pr for Au+Au 
collision at ^sjvjv=200 GeV in AMPT string melting model. 
The red dash curve is the fitting of data with the form of 
Eq. and the green dash-dot curve is a plot of Eq. (|10p . 
The blue curve is the flow-only simulation. 
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FIG. 4: (Color online)(a) (1 + Cpg) (solid points) and (1 + 
C'pb) (red dash curve) from AMPT are divided by (1 ~f C#b) 
from flow-only simulation, (c) (?-parameter as a function of 
Pt- 



the automatic decrease of Cfb (dominated at large pt)- 
After fitting, we could obtain the first part of Eq. ([S]): 



c% = -q{PT)\ 



(10) 



where the effect of flow multiplicity correlation between 
the forward and backward regions has been canceled. 
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In Fig. [31 we present the fitting of AMPT data by 
tlie red dash curve. We can see that this form can 
well describe the trend of the data, and we obtain that 
{il^r]^) ~ 0.013. For comparison, we fit the particle dis- 
tribution with Levy function, which can well describe the 
data, and generates flow particles with the same distri- 
bution. As discussed before, the flow particles will ran- 
domly fall into each side in this simulation. The result 
of this flow-only simulation is shown as the blue curve 
in Fig. [31 From the plot, we can see that the simulation 
stays zero at low pt as expected since there is neither jet 
nor correlation. While at intermediate px range, both of 
the AMPT data and the flow-only simulation decreases 
as PT increasing. We found that the AMPT result in- 
cluding jet contribution is only slightly lower than the 
flow-only simulation. Therefore, we argue that the de- 
creasing trend of the forward-backward correlation func- 
tion as shown in [loj may not only due to the jet contri- 
bution, but mostly due to the artificial decrease for the 
particle-rare case. 

To see the details clearly, the results from AMPT 
model are divided by that from the flow-only simulation. 
Here we define 
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The deviation from unity of this ratio refiects the dilu- 
tion of (NfNb) AMPT due to pairs containing jets. From 
the plot, we can see that the ratio begins to deviate from 
1 around 2 GeV/c, which indicates that the jet contribu- 
tion begins to join in this region. The difference between 
the points and the curve comes from flow multiplicity 
correlation. 

For qualitatively study, we assume that the square root 
of the ratio is proportional to the particle yield which is 
expected to work for pT > 2.2 GeV/c as shown in Fig. [2l 
and finally g can be extracted from (1 — v^). Fig. [H^b) 
shows that the g-parameter is about 10% at 2 GeV/c 
and reaches 30% at 3 GeV/c. Our method works best 
for central collisions in which away-side production of a 
jet is suppressed thus those particles can be considered 
as being "melted" into fiow particles. While for the non- 
central collision, the g quantity will be reduced due to the 
not quenched away-side jets. In that case, the correlation 
function measures the extra jets in one side. 



IV. THE CONNECTION TO PREVIOUS 



OBSERVABLE 



The forward-backward elliptic anisotropy correlation 
has been discussed in Ref (ol. lTo|. This correlation func- 
tion also subjects to the small number statistic discussed 
before. Here we redefine the elliptic correlation function: 
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(12) 



where vf'"^'' refers to the sum of the V2 of particles in 
forward (backward) region in each event. The {V2) we 
measured in experiments is constituted of two parts: flow 



and jet w^: 



{V2) = il-g){vi)+g{vi) 



(13) 



describe the collective behavior of hydro, and V2 is 
related to the parton energy loss in different directions. 
The relation between and pure jet (flow) V2 can then 
be written as: 
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A{r,^V^){l-g)^{v^) 



(14) 



The differences between this and Ref [T3 is coming from 
the —1 in the definition. One should pay attention that 
{ij^rj^} here is not the same as in the multiplicity correla- 
tion. For multiplicity correlation, rj^^^^ refers to the total 
multiplicity fiuctuation in each side. While for V2 correla- 
tion, rj^^^'^ have azimuthal dependence, i.e., the in-plane 
and out-of-plane multiplicity fiuctuation may cause the 
statistic fiuctuation of V2- To be clearer, the Eq. ([M]) can 
be written as 
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Therefore, for the hydro dominate case, i.e., 17 — > 0, 

Cp^g — 4-^^. The Cp^g describes the V2 fluctuation 

at low pt, and the large value of the Cp^g observed in 
Ref [l3] could be understood according to the magnitude 
of the flow fluctuation studied in Ref [l^ . For interme- 
diate pt; however, the decreasing trend of observed 
are mostly due to the rare fiow particles as discussed 
before. Therefore, it is interesting to check the behavior 
of Cp^g in the LHC energy where particle yield is much 
higher. 



CONCLUSION 



In summary, we have proposed a forward-backward 
multiplicity correlation function Cpg to study the jet 
contribution. We find that the Cpg is sensitive to the 
jet contribution, and its value will be enlarged by the 
forward-backward multiplicity correlation. In the inter- 
mediate Pt, We find that the Cpg is sensitive to small 
number statistic and will automatically decrease even 
there are no jets involved. This effect should also be 
present in previous study of Cp^g , and might be misun- 
derstood as the contribution of jets. 

We study the Cpg in Au+Au collisions at ^/snn ~ 
200 GeV with AMPT string melting model. When pT 
is lower than 1.4 GeV/c, Cpg is independent oi pt, and 
{rj^f]^) is about 0.01. At intermediate pt range, Cpg 
decreases with pt, and the result from AMPT model is 
lower than the flow-only simulation of the same yield. It 
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indicates that both of the the jet contribution and the 
procedure of excluding events with empty bins will cause 
the decrease of Cpg . The estimated jet fraction is much 
smaller than previous study. 

Finally, we discussed the connection of our study to the 
forward-backward elliptic anisotropy correlation func- 
tion. The large value of Cp^g is found due to V2 fluc- 
tuation, and the decrease of Cp^g observed previously is 
caused by both of the jet and the tail of flow particle 



distribution, with the latter dominates. 
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